Traveling wave fronts and the transition to saturation 
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We propose a general method to study the solutions to nonlinear QCD evolution equations, based 
on a deep analogy with the physics of traveling waves. In particular, we show that the transition to 
the saturation regime of high energy QCD is identical to the formation of the front of a traveling 
wave. Within this physical picture, we provide the expressions for the saturation scale and the gluon 
density profile as a function of the total rapidity and the transverse momentum. The application 
to the Balitsky-Kovchegov equation for both fixed and running coupling constants confirms the 
effectiveness of this method. 
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'{C)'. I. INTRODUCTION 

o : 

Phenomenology of high energy deep-inelastic scattering in the framework of QCD leads to the study of evolution 
equations for the gluon distribution function Af(k, Y) at high density 0, 0, 0, !H ■ ^ is the rapidity of the 
^ ■ evolved gluons and plays the role of the evolution variable, and k is their transverse momentum. The main feature 
of these evolution equations is that they contain nonlinear damping terms which make them difficult to solve by 
general methods. These damping terms reflect the effect of saturation, which is due to the recombination of partons 
densely packed together. The evolution at low density is well-understood and is described by a linear equation. The 
deep saturation regime can also be evaluated. However, the transition between these two regimes is still a challenge. 
£NJ ' Some interesting approaches for solving non- linear evolution equations have been proposed . Recently, the use of 
solutions of the linear evolution supplemented by "absorbing" boundary conditions has been proposed to determine 
the energy dependence of the saturation momentum [6| . 

The purpose of the present paper is to propose a new method for handling this problem. It is based on a deep 
connection that we have found between the formation of traveling wave fronts in nonlinear physics and the transition 
to saturation. 

. We expect our method to be very general and applicable for various types of equations. In the present paper, we 
are going to apply it in detail to the Balitsky-Kovchegov (BK) equation [j], which reads 



in 
m 
O 



d Y M= a X (~d L )JV-aAf 2 , (1) 



pi ' where L — log(A 2 /AQ C - D ). The QCD coupling a is either fixed or runs with k, in which case it is given by 

•i—i . 

^ ■ The function X is the well-known Balitsky-Fadin-Kuraev-Lipatov (BFKL) kernel |jj 

. d; x ( 7 ) = 2^(1) - - ^(i - 7) (3) 

and x is an integro-differential operator which may be defined with the help of the formal series expansion 

X (-d L ) = x(7o)l + x'(lo)(-d L - 7ol) + ix"(lo)(-d L - 7o l) 2 + hX (3) (lo)(-d L - 7ol) 3 + ■ ■ ■ (4) 

for some given 70 between and 1, i.e. for the principal branch of the function x- The BK equation is expected 
to capture the main features of saturation effects in high energy scattering processes like virtual photon-nucleus (or 
proton under certain conditions) total cross section. 
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In a previous paper and in the case of fixed QCD coupling, we have already shown |8( the intimate connection 
between the limiting speed of a wave front and geometric scaling, namely the dependence of J\f(k,Y) on a unique 
scaling variable k/Q s (Y), where Q s is the so-called saturation scale. Interestingly enough, geometric scaling has been 
observed in the HERA data on inclusive 7* — p scattering 0- One aspect of the present letter is to go beyond these 
asymptotic results and to show how the approach to geometric scaling is related to the mechanism of wave front 
formation. This requires a method to get information on the form of the wave front when rapidity increases. Another 
goal is to extend our analysis to both fixed and running QCD coupling. 

In Sec[ni we present the general method that we propose. In Sec lIIII we treat the fixed coupling problem and in 
Sec lIVI the running coupling one. For both cases we derive the expressions for the saturation scale Q S (Y) and for the 
gluon density J\f(k/Q S (Y),Y) in the vicinity of the geometric scaling region L ~ log(Q^/AQ CD ). The final SeclVlis 
devoted to a discussion and summary. 

II. PHASE VELOCITY, GROUP VELOCITY AND WAVE FRONT FORMATION 

Let us outline the method how to investigate traveling wave solutions to nonlinear partial derivative equations. 

Our analysis comes from the physics of linearly unstable front equations of diffusion type. They appear e.g. for 
Rayleigh instabilities in hydrodynamical flows, directed polymers in random media, chemical and bacterial growth 
models etc.jlflj]. The evolution equations have the property to admit traveling wave solutions. Their existence is due 
to nonlinear damping while some of their large time characteristic features (front velocity, front shape) are determined 
by the linearized equation about the unstable state. The front is said to be "pulled" by this unstable state which 
appears as an unstable fixed point of the equation. The main role of the nonlinearities is to cause saturation. 

Let us illustrate these features on the generic equation 



where ED is a linear operator acting on functions of x, f(u) contains the nonlinearities, with the requirement that 
D • u(x = const., t) — 0, and /(0) = 0, so that u = is a fixed point. This fixed point is furthermore required to be 
unstable with respect to perturbations, which has to be checked in each particular case. 

The simplest case is the Fisher and Kolmogorov-Petrovsky-Piscounov (KPP) equation ^lj corresponding to B = d 2 
and f(u) = u — u 2 , for which rigorous mathematical results are available 01 • Using a change of variables, we have 
shown 8] that the BK equation reduces to the KPP equation when its kernel is approximated by the first three 
terms of the expansion (@J) . 

More generally, on physical grounds 1 , we expect the methods developed to treat the KPP equation to be valid for 
any saturation equation in QCD. Indeed, the common feature of these equations is that their nonlinearities cause 
saturation of the parton densities, and the fixed point of low gluon density is clearly linearly unstable under rapidity 
evolution. This instability is physically related to the growth of parton densities with rapidity given by the BFKL 
equation 0. 

We shall study the form of the wave front at large times by working in its own reference frame which is defined by 
the change of variable x = x WF + vt, and in the neighbourhood of the front, namely a; WF <C vt. The velocity v of the 
wave front will be determined by consistency as the velocity of the slowest wave present in the initial wave packet. 

The solution to the linear part of equation J5j corresponds to a linear superposition of waves: 



where a; (7) is the Mellin transform of the linear kernel (d 2 + 1 for the KPP equation) and defines the dispersion 
relation of the linearized equation. In particular, each partial wave of wave number 7 has a phase velocity 



d t u(x,t) = B • u(x,t) + f(u) 



(5) 




(6) 
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At variance with the case of t he KPP equation, mathematical proofs are lacking for properties of general nonlinear equations. However, 
a large number of results 10] have been found which e.g. numerically confirm the existence and properties of traveling wave solutions 
to equations of the type ISl . 



whose expression is found by imposing that the exponential factor in Eq.© be time independent for v = Vu,(j). By 
contrast, the group velocity is defined by the saddle point j c of the exponential phase factor 



(ko 
d~j 



(8) 



Let us show that the actual velocity v of the front depends on the competition between the singularities of ^0(7) 
and the saddle point of the phase factor. 

In Eq.©, 1*0(7) is the initial condition in Mellin space, i.e. the Mellin transform of u(x,t — 0). The physically 
appropriate assumption is that 1*0(7) is a monotoneous function smoothly connecting 1 to as a; goes from —00 to 
+00, with asymptotic behavior u(x, t — 0) ~ e~ loX Then ^0(7) has singularities on the real axis starting at and 
extending to the left, and at 70 to the right. A simple case study is u(0,x) = 1 for x < and u(0,x) — e~ loX for 
x > 0, whose exact Mellin transform is 1*0(7) = 1/7 + V(7o — 7)- The contour C in JHJ runs parallel to the imaginary 
axis and crosses the real axis at a point between and 70. There are three relevant cases (see Fig.l): 

(i) 70 < 7 C : dominance of initial conditions. The integral © defining u is dominated by the pole at 70, thus 
u ~ exp{— 7o(x WF + vt) +w(7o)i} and the velocity of the front is determined by consistency to be v = w(7o)/7o, which 
ensures the time- independence in the reference frame of the front. 

(ii) 70 > 7 C : dominance of the saddle point. The integral in Eq.© is now ~ cxp{— 7 c (a; WF + vt) + Lu("f c )t} and thus 
the velocity v of the front is defined by 

w(7 c ) . w(7) 

v — v„ = = mm . (9) 

7c 77 

Indeed, in this case, the group velocity is identical to the minimum of the phase velocity 2 . In the previous case instead, 
the front velocity is v = w(7o)/7o > v g- 

(Hi) 70 = 7 C : critical case. Both the pole and the saddle point contribute, and v = v g . 

Now that we have determined the velocity of the wave front, reporting it in Eq.@ and performing the integral in 
the large t limit gives the asymptotic expression for the form of the front: 

t— kx> I e 7c WF 11 70 > 7c 

The typical wave front solutions are illustrated on Fig. 2. The upper plot displays a typical feature of the wave front 
development for the case 70 < 7 C - Here the front moves with the velocity v~ > v g and keeps the memory of the initial 
condition. The lower plot illustrates the case 70 > 7c- In this case, the wave front develops with the fixed group 
velocity v g . The forward part of the curve (after the knee) is the part of the front which keeps the memory of the 
initial condition, which is a typical feature of a "pulled" front [lfj. The middle plot is for the critical case 70 = 7c- 

In the following, we will focus on the case in which the initial condition is steep enough, so that the velocity of the 
front is the group velocity © and the form of the wave front is the second equation of (|10|l . 

It is interesting for our purpose to go beyond the asymptotics (IK)ll and study how the front reaches its asymptotic 
form as a function of time. For this sake, we shall use the following ansatz in the neighbourhood of the front: 

u(x,t) = t a G( XwF+c{t) )e-^ x ^ +c( - t » , (11) 

* large V t a J 

where c(t) is a sub-dominant function of time (with respect to v g t). Indeed, the front velocity receives some sub- 
asymptotic corrections w.r.t. the constant asymptotic velocity v g , contained in the derivative c(t) — dc(t)/dt. The 
factors t a describe a diffusion-like evolution which will be made explicit in each particular case. Hence the function 
t a G represents a subasymptotic correction to the form of the front UlUfl at large times. This ansatz was proven to be 
relevant for KPP-like equations 13]. Its key property is to be a solution of the linearized equation at large t in the 
neighbourhood of the front, for some a and c(t) to be determined. Note that since u(x, t) only depends on x WF + c(t) 
when the front has reached its asymptotic form of a pure traveling wave, G(z) is required to vanish like z for z — ► 0. 

The features of this ansatz are displayed on the lowest plot in Fig. 2. One sees that a steep initial condition 
characterized by the exponential decay slope 70 evolves, after a given time t, into a front characterized by the critical 
exponential decay slope j c . Its width is of order t a as we will show later on. 



2 Note that the relation v g =v v , in analogy with wave physics, has been written in Ref. 
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FIG. 1: The front velocity as a function of the steepness 70 of the initial condition at x — > 00 (u(x,t = 0) ~ e w:c ). Continuous 
line: the actual front velocity. Dashed line: the phase velocity for a wave with wave number 70. The abscissa points 7_, 7+, -y c 
correspond to initial conditions resp. below, larger than and equal to the critical value (resp. cases (i), (ii) and (Hi), see text). 



Let us first apply the formalism exposed above to the BK equation with fixed coupling a. The linear part of Eq.JQ) 
is solved by 



The rapidity Y is the evolution variable, and thus plays the role of the time t, while L — log(/c 2 /AQ CD ) plays the role 
of the spatial coordinate x. 

Let us discuss the relevant initial condition Noil) f° r QCD. The essential property for the derivation is that at 
large L (i.e. large k), J\f(k,Y = 0) ~ fc~ 2 : this corresponds to the color transparency property of gauge theories, and 
yields 70 = 1. For L small or negative, the dynamics is unknown since it belongs to the nonperturbative regime of 
QCD. However, interestingly enough, it does not play a direct role in the final result for the front formation. One may 
assume that J\f(k,Y = 0) is bounded |8§, which is a way to impose unitarity on the initial condition. In Mellin space, 
this leads to singularities of Afo(j) distributed on the intervals ] —00, 0] and [1, -+00 [. We shall verify in a moment that 
the wave number j c which minimizes the phase velocity, i. e. the wave number of the partial wave which moves with 
the group velocity, is between and 1. Thus we are in the case (ii) of Sec[n]in which the saddle point dominates. 

We shall now proceed to the determination of the analogous of the group velocity v g . By comparison of Eq. Ijl2|l to 
Eq.lJHJl, one sees that the dispersion relation reads 



III. BK EQUATION WITH FIXED QCD COUPLING 




(12) 



The phase velocity 1^(7) 



^(7) = axil) ■ 

w(7)/7 has a minimum at 7 C . From Eas. Q8l9(l . 7 C solves the implicit equation 

IcX'ilc) = Xilc) , 



(13) 



(14) 
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FIG. 2: Picture of the evolution of a front for three typical initial conditions. Dashed line: the initial profile u(x,t = 0). Full 
line: the evolved front after at large time t. Upper plot (case (if): the exponential slope 70 of the initial condition is 7_ < "/ c . 
The front propagates at the phase velocity V- = 1^(7-). Middle plot (case (Hi)): 70 = 7 C . The front propagates at velocity v g . 
Lower plot (case (ii)): 70 = 7+ > 7c The front still propagates at velocity v g , the size of the front is finite of order t a , and its 
slope is given by 7 C . 

which yields 7 C = 0.6275... . This number indeed stands between and 1, and thus we are in a similar situation as 
the one depicted on the third plot of Fig. 2: the front velocity is blocked at the value 



To study the formation of the front, we will stick to the "diffusive" approximation of the BK equation JIJ, defined 
by keeping terms up to the second order only in the expansion 



which is equivalent to the KPP equation 8]. 

We now use the ansatz as the appropriate solution of the linearized version of Ea. l|17() near the wave front. 
Setting u=M , x = L and t=Y, the linear part of Ea. 1)17(1 turns into an ordinary differential equation for G(z), namely 




(16) 



d Y Af = -v g d L N + ±a X "(lc){d L + 7d) 2 Af - aN 2 , 



(17) 



|S X "( 7c )y- Q G"(z) + (azV 



c(Y))G'(z) + Y a -\j c c{Y) Y - a)G{z) = 



(18) 



where z is related to the physical variables through 

+ \ (T.— n V^-i- r(V\ 

(19) 



x WF + c(t) (L- v g Y) + c(Y) 



t a Y a 

The different terms in Ea. (|18|) contribute equally if a = 5 and c(Y) = (3/Y, where (3 is a constant. This value of a 
ensures that the terms proportional to Y~ a and Y a_1 contribute to the same order in Y . Using these values, Ea. ll8|) 
boils down to a simple second order differential equation 

ia X "( 7c )G"(z) + izG'(z) + (/? 7c - \)G{z) = . (20) 

The constant j3 is determined by requiring that G behave smoothly for z — > +00, which fixes uniquely 3 (3 = 3/(27 c ) 
|l3j |. The solution which vanishes at 3 = is unique up to a constant and reads 



G(z) = const, x i / — - z exp f — - — — - I . (21) 



In terms of physical variables, the gluon density near the wave front thus reads 



with the saturation scale defined as 



4 



7c 27, 

Qo absorbs undetermined constants but remains of order Aqcj> 



2(y) = Qgexp(a^r- ^-logF) . (2:!) 



IV. BK EQUATION WITH RUNNING QCD COUPLING 

We turn to the case in which the coupling in Eq.QJ runs with L as given by Eq.©. The BK equation reads 

bL dyN = X (— <9l) A/" - N 2 , (24) 
with b defined in Eq. J2J) . Using a double Mellin transform, the linearized version of the equation gives the solution 

N{k > Y) = I ^Ll l^ 7 ' w) exp {~ iL +u;Y+ iL x{i) ) ' (25) 



where 



*(7)= fdixil') • (26) 



At this stage, 7 is an undetermined constant. The initial condition A/"o(7, u>) depends now both on u> and 7, but has 
to satisfy the same requirements as in the fixed coupling case (color transparency at large k, saturation at smaller k) 
resulting in a specific analytic structure in 7 as already discussed in Sec. III. 

Performing the integral over uj in the large Y limit, the integral is dominated by a saddle point at uj c = y/bX(^y)/Y. 
The solution reads 



M(L,Y)= J |7-AAo(7)exp ^-yL + W 



4* (7) 



(27) 



3 The full asymptotics of G in the framework of the KPP equation can be found in Ref.|l3j. 

4 Note that 1231 reproduces the result of |i| after expansion of 7 C around i . 



where 7Vo(7) keeps track of all prefactors. Now interpreting vF as the time variable t, the gluon density turns out to 
be also a linear superposition of traveling waves labelled by 7, and thus can be treated by the general method exposed 
in Sec[n In this case, the asymptotic dispersion relation is 



-(7) = ^ ■ < 28 > 

The minimum 7 of 1^,(7) = uj(^)/~i satisfies the equation 

2X(7)= 7X (7), (29) 
still depending at this stage on the choice of 7 in Eq.JJHJ). The velocity of the front (see Eq.JHJ) is given by 



v a = ]J-£-- (30) 

Requiring the velocity v g to become independent of the choice of 7 leads to the constraint dv g (j)/drf = = dv g (j)/drf, 
namely 

7X'(7)=X(7). (31) 

The comparison with Ea. ll4f) shows that 7 is identical to j c = 0.6275..., the critical wave number in the non-running 
case, and thus the velocity of the front reads 



v g = v(j c ) = y-^T^ ■ (32) 

Once the asymptotic velocity is known, we compute the transition to saturation through subleading corrections to 
v g . To this aim, one expands the linearized BK equation around 7 C to get 

b ^d t N = - b -^d L N+ [dlN + 2 lc d L N + 7 » ■ (33) 

One uses the ansatz for the solution to this equation at large t, which turns it into an ordinary second order 
differential equation. After tedious but straightforward compilation of all terms, one determines the coefficient a and 
the function c(t) by requiring that they all contribute equally. One finds a = § and c{t) = f3t~ 2 /' i , where (3 still needs 
to be fixed. The value of a comes from keeping relevant terms in Y~ a and Y 2a ~ l , which all contribute at large Y. 
Finally, Ea. (|33|) reduces to the Airy equation 

Q n {z) ^b^Vg_ 

X "(7c) 

There are two independent solutions to this equation, e.g. Ai and Bi. Imposing the regularity at z = +00 selects 



JcVgb 



1/3 



G(z) = const, x Ai ( ( ) (z - 4/3) ) . (35) 



The condition G{z) ~ z for small z is satisfied if 



, = .U^M)'" (1 , (36) 



4 V IcVgb 

where £1 = —2.338... is the rightmost zero of the Airy function. The complete result reads, in physical variables 

1/3 



where the saturation scale is given by 



(37) 




up to a multiplicative constant. 
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FIG. 3: Evolution of the reduced front profile in the case of fixed coupling (left) and running coupling (right). The reduced 
front profile (k 2 /Q 2 ) Jc JV(k/Q s (Y),Y) is plotted against log(k 2 /Q 2 ) for different rapidities. The various lines correspond to 
rapidities from 2 (lower curves, full line) up to 10 (upper curves). Note the similarity of the wave fronts, but the quicker time 
evolution (in y/t) for fixed coupling, by contrast with the slow time evolution (in i 1 ' 3 ) for the running coupling case. 



V. DISCUSSION AND CONCLUSION 



Let us discuss our final results for fixed 1221231) and running (|37I38|I coupling. In both cases, we get the leading 
exponential terms (resp. a(x(lc)/jc)Y and y/Hxilc) /^HcW) in the expressions for the saturation scale found in 
Refs. |6lll4| . These terms are supplemented by subleading corrections in rapidity, which are precisely those present 
in Ref.p, and which are respectively proportional to logF and Y" 1 / 6 . They are related to the slowing down of the 
velocity due to formation of the wave front. We expect this slowing down to be a very general feature of the transition 
to saturation. The expressions for the front profiles exhibit geometric scaling violations which are due to the time 
[i.e. rapidity) delay for the formation of the front. 

Comparing to the results which were obtained from the linear evolution equation supplemented by absorptive 
boundary conditions we have confirmed these results for the running coupling case 5 , up to a constant multiplicative 
factor in the definition of the saturation scale Q s . For the fixed coupling case, the same is true within the geometric 
scaling region k ~ Q S (Y). Our formula (|22[) provides a hint on the geometric scaling violations in this case also. 
A stimulating and fruitful challenge for future applications is to understand why our method and the one of Ref.0 
appear to be equivalent for the treated problems. It is interesting to note that numerical simulations of the fixed 
coupling case [l6j and of the full equation with running coupling [17| confirm the validity of the analytic estimates. 

To summarize, in this paper, we have proposed a general method, physically appealing, which allows to determine 
both the saturation scale and the transition of the gluon density towards saturation for the solution of the Balitsky- 
Kovchegov nonlinear equation. The method is based on a direct analogy with traveling wave solutions of nonlinear 
equations of the KPP type. The saturation scale is related to the velocity of the wave front, and the transition to 
saturation of the gluon distribution to the diffusive type formation of the front. In particular, we have studied directly 
the solution to the nonlinear equation, instead of assuming specific boundary conditions on the linear one. 

As first applications of this method, using a general ansatz |13| , we have derived in a rather simple way the physical 
solutions to the BK equation at large rapidity for both fixed and running coupling, see Eas. (|22l23|l and (|37I38|I . The 



5 Note that an extension to the (resummed) next- to-leading order BFKL kernel has been found in Ref.llf^. It will be interesting to see 
how this result can also be obtained by our method. 



approach to saturation (leading to geometric scaling violations) corresponds to ordinary diffusion (namely with a 
diffusion length going like y/i) in the fixed coupling case and to anomalous diffusion (in i 1 / 3 ) in the running coupling 
case. These contrasted diffusion patterns towards saturation are made explicit in the evolution of the reduced front 
profile (fc 2 /Qs) 7c Af(k/Q S (Y), Y) shown on Fig.3. 

We think that the method that we have presented is general enough to allow for desirable extensions to the full 
content of saturation equations. 
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